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Abstract We study the three-dimensional Anderson model 
of localization with anisotropic hopping, i.e., weakly coupled 
chains and weakly coupled planes. In our extensive numerical 
study we identify and characterize the metal-insulator transi- 
tion by means of the transfer-matrix method and energy level 
statistics. Using high accuracy data for large system sizes we 
estimate the critical exponent as u = 1.6 ± 0.3. This is in 
agreement with its value in the isotropic case and in other 
models of the orthogonal universality class. 

Previous studies of Anderson localization Q in three- 
dimensional (3D) disordered systems with anisotropic 
hopping using the transfer-matrix method (TMM) |g, 
g,^, multifractal analysis (MFA) ^ and energy-level 
statistics (ELS) g show that an MIT exists even for very 
strong anisotropy. In Refs. [0,||] , we studied critical prop- 
erties of this second-order phase transition with high ac- 
curacy. Here we shall demonstrate the significance of ir- 
relevant scaling exponents for an accurate determination 
of the critical disorder Wc and the critical exponent v. 
Previous highly accurate TMM studies for isotropic sys- 
tems of the orthogonal universality class reported v = 
1.54 ± 0.08 §, z^ = 1.58 ± 0.06 0, i' = 1.61 ± 0.07, and 
v = 1.54 ± 0.03 [Ol, whereas for anisotropic systems of 
weakly coupled planes h' — 1.3±0.1 and i/ = 1.3±0.3 was 
found H. We emphasize that this variation in theoreti- 
cal values has its counterpart in the experiments where a 
large variation of ly has been reported with values rang- 
ing from 0.5 [|| over 1.0 |l|, 1.3 0, up to 1.6 ||. 
Possibly this experimental "exponent puzzle" [|lj] is due 
to other effects such as electron-electron interaction [|l5[ 
or sample inhomogeneities |p^, |l6| , [l7|| . 

A further important aspect of anisotropic hopping 
besides the question of universality is the connection to 
experiments which use uniaxial stress, tuning disordered 
Si:P or Si:B systems across the MIT ||l|,|l|,[l4[|l|l. Ap- 
plying stress reduces the distance between the atomic 
orbitals, the electronic motion becomes alleviated, and 
the system changes from insulating to metallic. Thus, al- 
though the explicit dependence of hopping strength on 
stress is material specific and in general not known, it is 
reasonable to relate uniaxial stress in a disordered sys- 
tem to an anisotropic Anderson model with increased 
hopping between neighboring planes. 

We use the standard Anderson Hamiltonian ||l| 
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with orthonormal states \i) corresponding to electrons 
located at sites i = (x, y, z) of a regular cubic lattice with 



periodic boundary conditions. The potential energies e^ 
are independent random numbers drawn uniformly from 
[—14^/2, W/2]. The disorder strength W specifies the am- 
plitude of the fluctuations of the potential energy. The 
hopping integrals tij are non-zero only for nearest neigh- 
bors and depend on the spatial directions, thus i^ can 
either be tx, ty or tz. We study (i) weakly coupled planes 
with tx — ty — 1, iz = 1^7 and (ii) weakly cou- 
pled chains with tx = ty ^ \ — ^, t^ — 1 with hop- 
ping anisotropy 7 G [0,1]. For 7 = we recover the 
isotropic case, 7 = 1 corresponds to independent planes 
or chains. We note that uniaxial stress would be mod- 
eled by weakly coupled chains after renormalization of 
the hopping strengths such that the largest t is set to 1. 

The MIT in the Anderson model of localization is ex- 
pected to be a second-order phase transition ||l8|,|9|. It is 
characterized by a divergent correlation length ^00 (W) oc 
\W - Wc\~'^ @. To construct the correlation length of 
the infinite system ^00 from finite size data S,m P,|2Q, 
^,[2^, the one-parameter scaling hypothesis [g3| S,m — 
f{M/S,oo) is employed. One might determine v from fit- 
ting to ^00 obtained by a FSS procedure |Q. Better ac- 
curacy can be achieved by fitting directly to the ^m data 
(9[|l^,|ll[. We use fit functions [|o) which include two 
kinds of corrections to scaling: (i) nonlinearities of the 
disorder dependence of the scaling variable and (ii) an 
irrelevant scaling variable with exponent —y (cp. Fig. nl). 
For the nonlinear fit, we use the Levenberg-Marquardt 
method |n,pQ|. The input data S,m for the FSS procedure 
are either (a) reduced localization lengths Am obtained 
by TMM with 0.07% accuracy and system widths up to 
17 X 17 for, e.g., the case of weakly coupled planes with 
7 — 0.9 B; or (b) integrated A^ statistics obtained from 
highly accurate ELS data (0.2% to 0.4%) and system 
sizes up to 50^ . 

When applying the TMM to our anisotropic systems, 
one has to consider two non-equivalent orientations of 
the axis of the quasi- ID bar: parallel and perpendicu- 
lar to the planes or chains. The localization lengths in 
the perpendicular direction are smaller than in the par- 
allel direction by a factor of about 1 — 7 for coupled 
planes and (1 — 7)^ for chains Q. The critical disorder 
Wc should not depend on the orientation of the bar [g[ . 
For strong anisotropics 7 > 0.9 this is difficult to verify 
numerically due to strong finite size effects as shown in 
Fig. |l|. By computing data for very large system sizes 
up to Af^ = 22^ (46^) for the case of weakly coupled 
planes with 7 = 0.9 (0.96) we can show that this finite 
size effect can be sucessfuUy modelled (cp. Fig. m) by an 
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Fig. 1 Am for coupled planes with 7 = 0.9 (per- 
pendicular orientation) with relative error 0.1%, W — 
7, 7.1, 7.2, • • • , 9.2 and M^ = 5^ 6^ 7^ • • • , 22^ The solid 
lines in the right part are fits to the data with y = 2.05±0.08. 
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Fig. 2 Results for Wc and v, for coupled planes with 
7 = 0.9, obtained from FSS of (parallel-direction) TMM data 
(open symbols) and ELS data. The error bars show the 95% 
confidence intervals. The dotted (dashed) lines represent the 
error bounds for u = 1.62 ±0.07 (1.45 ±0.2) of TMM (ELS). 
The solid line marks the result of llQj. The goodness of a fit 
is reflected in the size of the symbol. The 2 thick error bars 
mark high quality ELS fits for large system sizes. 

irrelevant scaling exponent and Wc is indeed the same 
for both orientations. 

In Fig. H, we show fitted values obtained by FSS 
of TMM data for different choices of expansion coeffi- 
cients in the nonlinear fit procedure. We conclude v = 
1.62 ±0.07 and Wc = 8.63 ±0.02. In Fig. |, we also show 
the results for FSS of highly accurate ELS data (0.2% 
to 0.4%) and system sizes up to N^ = 50'^. The error 
estimate is larger and the values of Wc and v are much 
more scattered than before. Comparing the spreading of 
the Wr and v values with their confidence intervals, the 



error estimates appear to be too small. E.g., the 95% 
confidence intervals of the smallest and largest Wc value 
do not overlap. We therefore estimate i/ = 1.45 ±0.2 and 
I^c = 8.58 ± 0.06 §. 

In conclusion, our results confirm the existence of 
an MIT for anisotropy 7 < 1 for weakly coupled planes 
found previsouly in studies using TMM |^, MFA Q, and 
recently by ELS H. We have shown that large system 
sizes, high accuracies 0,g] and irrelevant scaling expo- 
nents are necessary to determine the critical behavior 
reliably. Our results are in good agreement with other 
high accuracy TMM studies for the orthogonal univer- 
sality class |^,^|l^,|2J]. These numerical estimates seem 
to converge towards j^ « 1.6. 
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